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Abstract. The interaction between supernova ejecta and circumstellar matter, arising from previous episodes of 
mass loss, provides us with a means with which to constrain the progenitors of supernovae. Radio observations of 
a number of supernovae show quasi-periodic deviations from a strict power-law decline at late times. Although 
several possibilities have been put forward to explain these modulations, no single explanation has proven to be 
entirely satisfactory. Here we suggest that Luminous Blue Variables undergoing S-Doradus type variations give 
rise to enhanced phases of mass loss which are imprinted on the immediate environment of the exploding star as a 
series of density enhancements. The variations in mass loss arise from changes in the ionization balance of Fe, the 
dominant ion that drives the wind. With this idea, we find that both the recurrence timescale of the variability, as 
well as the amplitude of the modulations are in line with the observations. Our scenario thus provides a natural, 
single-star explanation for the observed behaviour that is, in fact, expected on theoretical grounds. 
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1. Introduction 

Massive stars and core-collapse supernovae (SNe) play an 
important role in driving the chemical evolution of galax- 
ies and in shaping the dynamics of the interstellar medium 
over all cosmological epochs, since the very first stars came 
into existence. Despite their importance, the lives and 
deaths of massive stars are poorly understood; in partic- 
ular, it is not known with any degree of certainty which 
massive stars produce which SNe. 

While progress is gradually being made in the direct 
identification of core-collapse SN progenit ors by search- 
ing for these in pre -explosion images (e.g. |Smartt] [2002: 
Ivan Dvk et al . 2003), such studies are limited to either ex- 
tremely nearby SNe, and/or those occurring in uncrowded 
regions. Perhaps as a result of the stellar initial mass func- 
tion, most direct detections of progenitors have yielded 
masses not significantly greater than ^IOMq. 

Although the evolution of the most massive stars (with 
M>30 Mq) is largely unconstrained, it is generally ac- 
cepted that mass loss drives these objects through the 
O star, Luminous B lue Variable (LBV), and Wolf-Rayet 
(WR) phases (e.g. IChiosi fc Maedeil llQSdh Significant 
progress has been made in recent years in our un- 
derstanding of the mass loss properties over the vari- 
ous evolutionary phases of massive stars. In a series of 
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papers, IVink et alJ l|l999t 1200(1 l200lh provided predic- 
tions of mass-loss rates of O and B supergiants, LBVs 
l)Vink fc de Koterlr2002lh an d late-type WR-type st ars as 
a function of metal content llVink fc de Koteiil200filh 

SNe that result from the core-collapse of massive stars 
explode in environments that have been modified by mass 
loss from the progenitor. The SN ejecta interact first with 
this circumstellar material before interacting with inter- 
stellar material. We might therefore expect that the dis- 
tinct differences in wind properties over the lifespan of a 
massive star might be imprinted onto the resulting circum- 
stellar media (CSM). Furthermore, we would then expect 
these differences to be manifested in the interaction be- 
tween the SN ejecta and surrounding material. By quan- 
tifying these differences one can then constrain the evolu- 
tionary phase of the exploding object. 

Over the last decades, radio observations of SNe have 
provided a means with which to constrain the density of 
the CSM around core-collapse SNe. The inferred mass- 
loss rates from modelling of most radio SN light curves 
yi eld values of M ~1 0~ 6 - lO _4 M0yr _1 (e.g. compilation 
m IWeiler et al.l l2002h Unfortunately, these average mass- 
loss rates are generally only accurate to within a factor of 
~10, and are typical of almost all types of massive star, 
making it difficult to pin down the evolutionary phase 
during which core-collapse occurred. 
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A small subset of radio SNe, however, show quasi- 
periodic modulations in their radio lightcurves. In what 
follows, we argue that this type of modulation may be 
the result of an LBV star that underwent S Dor varia- 
tions which entailed opacity changes in the wind-driving 
region, resulting in varying mass-loss rates. 

2. LBV mass-loss variability 

LBVs are unstable massive stars, located in the upper 
portio n (L/Lp> <; 5.5) of the Her tzsprung-Russell dia- 
gram IjHumphrevs fc David son 1994). They exhibit differ- 
ent types of variability which can be divided into three 
categories: (i) small-amplitude (~0.1 mag) "micro" vari- 
ability, which is common amongst blue supergiants, (ii) 
"moderate" S Doradus variations of ~1 - 2 mag (SD 
phases), and (iii) truly "giant" eruptions, of which the 
Galactic cases P Cyg and Eta Car are mo st well-known- 
Accor ding to an extensive list compiled bv lvan Genderenl 
()200l[) . the SD-phases occur on two timescales: < 10 yrs. 
("short SD" phases), and > 20 yrs. ("long SD" phases). 

Spectroscopically-determined LBV mass-l oss rates are 
found to ch a nge on S Dor tim e- scales (e.g. IStahl et all 
l200lL 1200.1) . IVink fc de Koterl ll2002l) attributed this 
mass-loss variability to the ionization and recombina- 
tion of the Fc lines which are responsible for driving the 
wind. This leads to "bi - stable" behaviour in the wind 
dPauldrach fc Puldll990t lLamers et aflll995l IVink et all 
1999) causing the star to flip back and forth between two 
states: that of low mass loss, high- velocity winds, to high 
mass-loss, low velocity, winds. This bi-stable behaviour re- 
sults i n jumps in the mass-los s rate of a factor of ~2 for 
LBVs l|Vink fc de Koterl l2002). and a decrease by a fac- 
tor of ~2 in the terminal wind-v elocity between O and B 
supergiants ijLamers et alJll9*9*^) . The wind density would 
therefore be expected to change by a factor of ~4 on the 
timescale of the S Doradus variations. In the absence of 
any other material around the star, this would result in a 
pattern of concentric shells of varying density. 



3. Radio SN Lightcurves and stellar wind 
interaction 

The properties of Radio SNe (RSNe) lightcurves and the 
model for the interaction of the SN blastwave with the 
surrounding CSM as described by JCheyalierl l|l982albl ) 
has been reviewed by IWeiler et al.1 l)2002f) . In short, the 
source switches on, reaches a peak, and generally decays 
smoothly over a timescale of years. The peak is attained 
at different times (days to months) for different radio fre- 
quencies, depending on the optical depth of either external 
and/or internal absorbing material. For a discussion on the 
relative importance of free-free absorption (FFA) by an 
external - possibly clumpy - medium, versus sy nchrotron 
self-ab sorp tion (SSA), we refer the rea der to IChevalierl 
l|l998l) and lFransson fc Bidrnssonl i 1998). The radio emis- 
sion nevertheless places constraints on the density p, and 
therefore on the ratio of M to the terminal wind velocity 



(u w ) of the CSM of the SN progenitor: p oc M/v w r 2 . Using 
RSNe lightcurves, and generally assuming a typical i> w for 
a red (super)giant (~10-20km s _1 ) constraints on M of 
pre-SN wind s have been placed on about a dozen RSNe 
(see Table 3. IWeiler et alJl2002|) with typical M values in 
the range lO^-lO" 4 il^yr -1 . 

The shock radius can be constrained through radio ob- 
servations, whilst the measured variability timescale is set 
by the wind velocity of the progenitor (see below). Ideally, 
one would wish to spatially resolve the shock radius di- 
rectly, as was possible for SN 19 79C using VLBI tech- 
niques jBartel fc Bietenhol d E)03) . Usually, one needs to 
assess the relative importance of SSA vs. FFA. If SSA pre- 
vails, and if the peak flux is available, i.e., at the transi- 
tion between the optically thick to thin phases, assuming 
equipartition between the magnetic fields and relativis- 
tic electrons yie lds an estimate of the size of th e radio- 
emitting region l|Shklovskiilll985t ICheva licr 1998). While 
if FFA is dominant, it is usual to simply adopt the mea- 
sured velocity from optical lines at early times together 
with parameter s from modelling th e light curves at differ- 
ent frequencies I Weiler et al.ll2002^) . 



4. Modulations in radio SNe lightcurves 

Variations in the radio lightcurves have been reported 
for about a dozen SNe of both type II and type Ibc 
(see Table 4 of Soderberg et al. 2005). These can gen- 
erally be interpreted as due to variations in the mass- 
loss rate due to a changing evolutionary phase of the SN 
progenitor. A subset of radio SNe show epi sodic bumps. 
Thes e include the class ic cases of SN 1979C l|Weiler et alJ 
Il992l) and SN 1998bw feulkarni et al.lll998l) tha t was as- 
sociated with GRB 980425 l|Galama et all 119981). as wel l 
as two recent r adio SNe, SN 2001ig llRvder et al.l l2004|) 
and SN 2003bg llSoderberg et al.ll2005l) . This has brought 
the number of radio SNe that show recurrent modula- 
tions in their radio light curves, to four. A caveat to the 
above is that finer sampling of radio light curves may re- 
veal the extent to which such behaviour is commonplace. 
Interestingly, SNe 200 lig and 2003bg are strikingly simi- 
lar both in terms of the amplitudes and timescales of the 
variations (see Fig. 8 in Soderberg et al. 2005 for a com- 
parison). We focus our discussion on these two SNe, and 
defer a discussion of SN 1979C and 1998bw to 



4.1. SN 2001 ig 

SN 2001ig is a transition object: on the basis of opti- 
cal spectroscopy, it was initially classified as type II (i.e. 
showing H lines) bv lMatheson fc Jhal ( 200ll) but metamor- 
ph osed into type Ib/c object (i.e. no H lines, weak He lines) 
bv lFnippenkof(|2002j) about 9 months later. This suggests 
that it has lost most of its H-rich envelope l|Rvder et alJ 
120041) . We will ret urn to this po i nt in jCQ 

For SN 2001ig iRvder et all l|2004l) used the model of 
IWeiler et al.l |2002) to derive parameter fits to the radio 
lightcurves yielding parameters such as the radio spectral 
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index a, the decline rate of the optically thin phase (3, 
and importantly, at least one parameter 5 (= —2.56) that 
describes how the optical depth of the local CSM changes 
with time. The relevance of 5 is that it relates to the rate at 
which the SN blast-wave radius rew changes with time: 
trw oc t m with m = —5/3. In the mini-shell model by 
IChevalierl (Il982albh . m = 1 implies a non-decelerating 
CSM. For SN 2001ig, m = 0.85. 

The recurrence timescalc (i) of the bumps in the radio 
flux of SN 2001ig i s deriv ed to be 150 days. Using Eq. (13) 
from lWeiler et all l|l98rt : 



AP 



shell 



dejecta U ( t 
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where ti is set by when the ejecta velocity was mea- 
sured. For SN 2001ig U -1 4d and w C j CCta ~ 15000 km s" 1 
l|Clocchiatti fc Prietol EoOl^) . We stress that the depen- 
dence on ti is weak. Assuming v w = 1 0-20 km s~ 1 . t ypica l 
wind velocities for red (super)giants, iRvder et al.l {2004) 
found a period P between successive mass-loss episodes 
that was too long for red (super)giant pulsations (typically 
a few 100 days), and too short for thermal pulses (10 2 -10 3 
years). For these reasons, they invoked an edge-on, eccen- 
tric binary scenario, consisting of a WR-star and a massive 
companion. The configuration of the binary system would 
lead to a pile-up of mass, giving rise to modulations in the 
radio light curve. Such a scenario was proposed on the ba- 
sis of hydrodynamical simulations bv lSchwarz fc Pringld 
l|l996h to explain the radio lightcurve of SN 1979C. 

One of the main differences between LBV and red gi- 
ant winds is tha t LBV winds ar e faster, typically 100- 
500 kms" 1 (e.g. iLeithererl Il997h c.f. 10-20km s" 1 for 
red (super)giants. This implies that if the progenitor of 
SN 2001ig were an LBV, the expected period between 
successive mass-loss episodes would be AP — 25 yr (for 
an assumed u w =200km s _1 ), consistent with the long 
SD phase. Inconsi stencies in the derivation of AP in 
IRvder et alJ l|2004j) results in their inferred radial spac- 
ings between the shells being too small (~ 2 x 10 15 cm, 
with ti = Id c.f. 1.6 x 10 16 cm from above) and might 
therefore be incompatible with the kind of WR pinwheel 
nebula they describe. 

4.2. SN 2003bg 

SN 2003bg was class ified as a type Ic (no H, He, Si lines) 
by iFilippenkol ((2003) ; in less than a month, the spec- 
trum of SN 2003bg had evolved to a type II SN show - 
ing broad P Cygni line s of hydrogen l|Hamuv et al Jl200^) . 
ISoderberg et alJ l)2005|) are able to fit the the radio light 
curve under the assumption of SSA. 

Remarkably, the short-term variability in the multi- 
frequency radio lightcurves of SN 2003b g, is found to be 
on ver y similar timescales as SN 2001ig. ISoderberg et alJ 
( 200^) infer enhancements of a factor of —2 in density dur- 
ing the deviations from pure power-law evolution. They 
consider a range of possibilities that might account for 
the observed modulations in the lightcurve, but favour a 



single-star progenitor model of a WR star that underwent 
episodes of intensified mass loss. However, they do not 
specify the nature of the physical mechanism that might 
give rise to such periods of enhanced mass loss. Using the 
measured expansion velocity for SN 2003bg at a time i; 
-20 days (|Filippenkoll2003t lHamuv et al.ll2003|) . we find 
AP — 25 years, again, consistent with the long SD phase. 

4.3. The cases of SN 1979C and SN 1998bw 

As outlined above, the relevance of SNe 2001ig and 
2003bg-like events is that there is evidence for regularly 
changing mass loss of the SN progenitor. These observa- 
tions may finally open up the possibility of yielding more 
stringent constraints on the pre-SN mass-loss properties 
than available from previous radio data. To our knowl- 
edge, the only massive stars known to fulfil the criterion 
of quasi-regularly spaced mass-loss episodes are LBVs. 

The only other SNe that have shown episodic mod- 
ulations in their radio li ghtcurves are SN 1979C and 
SN 1998bw. For SN 1979C. IWeiler et alJ l|l992h found a ra- 
dio modulation on a timescale of —1575 days, i.e. a factor 
—10 longer than for SNe 2001ig/2003bg. However, when 
we account for the rate at whi ch the blast-wave ra dius 
changed with time, 5 = -2.94 ijMontes et al.l l2000l) . we 
find A P — 160 yrs for Dejecta ~ 8000 km s~ 1 l|Panagia et all 
Il980l) . Now the circumburst radius of SN 1979C has been 
resolved with VLBI. Using the value thus de rived of the 
radio size (Fig. 5. iBartel fc BietenholdEoblj) . combined 
with the epoch of the radio modulations at —4 and 8 
years, we find a difference in an gular size of —0.5 mas . 
Assuming a distance of 15.2 Mpc ( Freedman et a,lJl200lfl . 
then implies mass loss variability on timescales of —150 
years (v w — 200 km s _1 ), i.e. consistent with AP derived 
above, but lon ger than for SNe 2001 ig/2003bg. We note as 
an aside that Ivan Genderenl (|200l[) does list some LBVs 
with SD phases that are greater than —100 years. 

Interestingly, for the more exotic case of SN 1998bw, 
a sim ilar analysis as in S I4.1I with m = 0.78 l|Weiler et alJ 
120011) . and a timescale of t ~ 50 days, we find 
AP — 6 years, which is consistent with the short SD 
timescale. We note however that the episodic bumps in 
the lightcurve of SN 1998bw are not as well-defined as the 
sinusoidal variations of SNe 2001ig/2003bg. 



5. Discussion 

We have proposed that LBVs undergoing S Doradus- 
type variations are viable progenitors of SNe that show 
quasi-sinusoidal variations in their radio lightcurves. The 
main motivation for this suggestion is that S Doradus- 
type winds provide a natural and physical mechanism that 
could give rise to the observed radio modulations. Both 
the timescales and the amplitude of the variations in CSM 
density as a result of the S Doradus variations are in line 
with the constraints imposed by the radio SN data. Our 
proposed scenario therefore provides support for a single- 
star progenitor system. Admittedly, there are uncertain- 
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ties associated with current analyses of the timescales of 
the radio SNe described in this paper, as well as with 
LBV S Dor variations in general. Until a significant body 
of radio SNe with episodic modulations becomes available, 
statistical considerations are deemed premature. If shorter 
SD variations are superimposed on the longer ones, then 
these might be revealed by further increasing the cadence 
of radio observations to a nightly basis. 

It may be relevant that both SNe 200 lig and 2003bg 
are transitional objects. SN 2001ig is an example of the 
group of lib SNe, that switch from type II (H and He 
both present) to Type lb (He present, H absent), whilst 
2003bg was first classified as a type Ic, but then underwent 
a change to a type II. These facts suggest that SN2001ig 
and 2003bg have a rather limited amount of H. This is also 
consistent with an LBV scenario, as LBV atmospheres are 
H-rich compared to WR stars, but He-rich compared to 
OB stars and red supergiants. 

We do not exclude the possibility that objects like SNe 
2003bg and 1998bw are in fact "stripped core" WR stars 
by the time they explode. However, given that WR stars 
are not known to exhibit the type of variability that we 
consider here, we believe it is more likely that the observed 
quasi-sinusoidal variations in their radio lightcurves are 
attributable to their prior LBV phase. This would mean 
that these objects are WR stars for only a very brief 
amount of time before they explode, and to all intents 
and purposes, we can categorize these objects as LBVs. 

The fact that our proposed scenario is simple and 
workable doe s not mean other s cenarios need necessarily 
be excluded. iRvder et al.1 l|2006h recently found evidence 
for a binary companion of SN 200 lig. This in itself is not 
surprising; it is well-known that a lar ge fraction of mas - 
sive stars occur in binary systems Ce.g. lMason et al . 1998). 
The more relevant question of whether the potential bi- 
nary companion ever had any direct influence on the star 
that exploded as SN 2001ig is more difficult to answer. 

The binary scenario th at could give rise to a pin- 
wheel nebula as proposed bv lRvder et alJ l|2004ll200^ and 
ISchwarz fe Pringiel l)l996|) requires the binary to be eccen- 
tric and viewed close to edge-on. Given the striking simi- 
larity of SN 2003bg to SN 2001ig. lSoderberg et alJl20o3l 
argue that the Wolf-Rayet pinwheel nebula hypothesis is 
unlikely, as the configurations of both systems must neces- 
sarily be very similar. They therefore propose a single WR 
progenitor, but do not put forward a mechanism for the 
CSM variability. The LBV progenitor scenario presented 
here may remove this shortcoming. 
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